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PREFACE 


The work described in this report was performed by the Propulsion 
Division of the Jet Propulsion Laboratory. 
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ABSTRACT 


This report describes the results of an initial effort at Jet Propulsion 
Laboratory (JPL) to nno lel the control loops of a 30-cm diameter electron 
bombardment thruster ? nd a transistorized power processor so as to predict 
its operation. Data, from which the model is made, are presented as well 
as comparisons between the computer outputs and test data from the JPL 
Solar Electric Propulsion systems laboratory. 
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I. INTRODUCTION 

Ion thrusters are at the state of development where they are being 
considered for use for primary spacecraft propulsion (Refs. 1-3) and for 
earth-orbital applications (Refs. 4 and 5). The effort at JPL has been 
directed toward the use of 30-cm diameter ion thrusters for primary pro- 
pulsion on a solar electric spacecraft. A computer simulation of the thruster 
control system, therefore, became of increasing interest in order to: 

(1) examine thruster operation as a function of expected variations in control 
loop parameters and thruster environment, (2) construct a rrxodel of the 
thrust subsystem for use in mission operations, and (3) examine interactions 
between the thruster and the solar array. 

The control of ion thrusters can, in general, be separated into two 
categories. The first deals primarily with the control of the thrust level 
after steady-state operating temperatures are reached. The second deals 
with bringing the thruster on line and monitoring its operation. The first 
category is covered in this report. The second, which uses a digital com- 
puter for timing and event monitoring, is discussed in Refs. 6 and 7. 

The work described in this paper represents the initial stages of the 
effort to simulate the thruster control loops. The m.odel was constructed 
using characteristic data, as available, from the latest 30-cm ion thrusters 
supplied from Hughes Research Laboratories for testing at JPL. Data was 
obtained initially from a modified 400 series thruster (Ref. 8). Some of 
the thruster characterizations were later updated to include data from a 
700 series engineering model thruster (Ref. 9). These data were used 
along with the control functions incorporated in a transistorized power pro- 
cessor of the type described in Ref. 10 to construct a model of the operation 
of the thruster control system. These loops were analyzed to determine the 
extent of the effect of control- loop parameters and the thruster thermal 
environment on control stability and accuracy. The model an'ilysis places 
the most emphasis on the very low frequency effects introduced by *he 
thruster. It does not include the high frequency oscillations normally present 
in the thruster discharge, since these oscillations are beyond the control- 
loop response frequencies. 
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II. THRUSTER OPERATION 


Figures 1 and 2 are cutaway views of the thruster modeled. I iquid 
mercury propellant is vaporized in separate feed systems for the propellant 
distributor, ion source chamber cathode, anc neutralize r. At present, 
thrusters operate with about 88% of the flow supplied to the propellant dis- 
tribution manifold, 10% to the discharge cathode, and 2% to the discharge • ' 

neutralizer. Electrons are drawn out of the high-density mercury plasma 
formed in the discharge cathode by the keeper electrode. These electrons 
are then dispersed by the baffle plate and accelerated into the main discharge 
chamber by a positive anode voltage. When in the main discharge chamber, 
these energetic electrons ionize the neutral propellant mercury and form a 
m rcury plasma. The divergent magnetic field, created by the magnets and 
soft-iron pole pieces, traps the ionizing electrons and enhances the ionization 
process Ions formed in the discharge chamber drift toward the screen elec- 
trode. '7hen passing through the plasma sheath, the ions are accelerated 
through the concentr'c holes In the screen and accelerator electrodes by the 
applied electric field. The ion beam ib then decelerated by the space charge 
forces in the ion beam to a potential slightly higher than the ambi.ent-space 
plasma potential and neutralized by the electrons emitted by the discharge 
neutralizer. Over a relatively wide range of thruster performance, thrust 
is directly proportirnal to the mass flow rate of vaporized propellant into 
the thruster. The thrust, T, is a function of ion beam current and screen 
voltage as shown by the equation: 

T = KI„ fV ~ 

B ^ s p 

where 

I- = ion beam current 
rJ 

K = constant 

= screen voltage 
Vp = plasma voltage 
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To throttle the ion engine, the flow rate of propellant into the thruster can 
be varied or the screen voltage can be varied. Propellant throttling is the 
normal mode used in order to simplify the power processor design. 

The neutralizer is quite similar to the cathode used to supply elec- 
trons to the main discharge chamber; however, it requires a much lower 
flow of mercury vapor to operate. 

The thruster is surrounded by a cover at the spacecraft-frame 
potential to prevent" the space plasma or neutralizing electrons from being 
drawn to the high-positive potential of the thruster body and creatirn 
severe power losses. 

There are three control loops required to operate the thruster. The 
main vaporizer is used as the primary control on the thrust level. Beam 
current is sensed and used to control the power to heat the main vaporizer. 
Th«* discharge power is also controlled. Control of the propellant- 
utilization factor (the ratio of the mass of the ions expelled to the propellant 
input mass) is accomplished by setting the discharge current level and the 
discharge voltage is controlled to a fixed level to limit the ion erosion rate 
at the cathode. Thirdly, the neutralize ’ ion beam couples to 'lie ambient- 
space plasma through a potential drop of several volts. This potential drop 
is indirectly controlled by regulating the keeper voltage by adjusting the 
neutralizer flow rate. 

A magnetic baffle is used to affect the impedance between the cathode 
and anode, thereby influencing the discharge stability. This magnetic baffle 
around the cathode is clearly shown in Figs. 1 and Z. The magnetic field 
from this coil is used to restrict the flow of electrons from th^* cathode to 
the anode and thereby alter the main discharge chamber direct current 
characteristics (this eff< ct is no+^ built into -^he present model because of the 
limited amount of data available at this time). 

One other point of .nterest can be seen in F'igs. 1 and Z, The main 
vaporizer drives an isolator followed by a propellant distributor, while the 
cathode vaporizer drives an isolator anrl tubing *o thi- centrally located 
cathode. All of these volumes introduce transport lags into the loops and 
tend to introduce low-frequency oscillations. 
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Figure 3 is a block diagram of the proposed Power Processor (PP) 
required to operate the thruster and the interconnections within the thruster. 
The primary purposes of the PP are: (1) to interface between a high- 

impedance sources and the thruster, (2) to provide all power supplies required 
by the thruster, and (3) to protect both itself and the thruster froi i electrical 
failure. There are twelve supplies required to operate and start the thrus‘.;er 
(not all on at the same time). The screen supply is the major supply 
(2 A at 1100 V) with the discharge supply next (12 A at 37 V). All of the 
rest (except for the accelerator supply) are low-voltage, low-power supplies. 
The accelerator normally runs at 500 V, 0. 004 A. Normal operating condi- 
-.ions for the thruster are listed in Table 1. 
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III. CONTROL MODEL 


A simulation of the control loops for a 30-cm diameter ion thruster 
is being developed at Jet Propulsion Laboratory in the form of a computer 
model. The simulation uses Continuous System Simulation Language (CSSL)\ 
and closes the loop by iteration between sets of measured data. The thruster 
and power processor supply characteristics are included within the model. 

Figure 4 shows a block diagram of the control components in a thruster- 
power processor set. This diagram is broken up and presented in a greater 
detail in Fig. 5. Most empirical thruster data used in the moJel were 
derived using a modified Hughes S/N 40< thruster. Some of the 400 series 
data has been updated by using information obtained with a Hughes S/N 703 
thruster. A description of these thrusters is contained in Table 2. The 
characteristics of a transistorized power processor of Ref. 10 are used. 

This unit was modified for operation with the 30-cm thruster. 

The model contains three major luops. These loops are: (1) the ion 

beam current-main vaporizer loop which controls primarily the thrust level, 
(2) the discharge voltage-cathode vaporizer loop which controls the discharge 
chamber voltage, and (3) the neutralizer keeper voltage-neutralizer vapor- 
izer loop which controls both the amount of mercury flow through the 
neutralize ind the potential o the exhaust plume. The first two loops (the 
beam currt ..t and the dischax'ge voltage) interact through the mass flow that 
each provides to the discharge chamber and by the power supplied to the 
discharge. All vaporizers are affected by the thermal environment and this 
effect is included in the model. The effect of variations in the magnetic 
baffle current, screen, voltage, and acceleration \'oltage is not included at 
this time. The computer-program listing for the model is inc’uded in 
Appendix A. Additional details of each of the 'hree loops are described in 
the following sections. 


programming Sciences Corporation, l.os Angeles, CA. 
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A. 


ION BEAM CURRENT -- MAIN VAPORIZER LOOP 


This loop is shown in Fig. 5a. The beam current is determined by 
the total mass flow rate and the power supplied to the discharge. This loop 
controls ihe main vaporizer flow rate and is affected by the cathode vapor- 
izer flow rate from the discharge loop. The empiiical relationship for the 
thruster characteristics as used in this loop is shown in F’^. 6, where the 
total mass flow- rate is expressed as the sum of tha main and the cathode flow- 
rates (rh-, f m^). A reference level for the thrust level is supplied to this 
loop in the form of an ion beam reference (1^ ref. ), The error signal is 
acted on by the vaporizer heater controller to heat the main vaporizer. This 
controller, whose characteristics are shown in Fig. 7, is normally a propor- 
tional devdce that saturates at large error levels. This controller deter- 
mines the value of the main vaporizer heater current. The applied heater 
power controls the vaporizer temperature. The vaporizer has a thermal 
transfer function G^(S) which can be represented by a first-order lag with a 
time constant of 250 s. 

The resulting vaporizer temperai-ure rise, /hen added to the tempera- 
ture of the vaporizer environment (TE 2), determines the vaporizer opera- 
ting temperature. Experimental results tor determining the vaporizer lip.ie 
constant are presented in Fig. 8. The heating and cooling time constants 
are very close as can be seen by the slopes of the normalized heating and 
cooling curves. Typical vaporizer curves arc shown in Fig. 5. The 
relationship between the vaporizer heater power and the rise ir /aporizer 
temperature ij presented in Fig. 10. 

Vaporized mercury lea\dng the porous tungsten surface flows into the 

discharge chamber througi. a high-voltage isolator, and then through a 

manifold which introduces some delay time. This main 'aporizer flow rate 

affe«^ts ihe discharge loop by (hanging the effective cathode mass flow rale, 

N1 , as follows: M = m t Kin., where K is an experimental) v determined 

e’ c c M ^ 

weighting factor. 

The ion beam (.uirent, as can be seen in P'ig. (>, is nc^f a monotonic 
function of total mass flcnv rate. The control loops will not function properly 
when th(= slope of tlie c u rv-e is negativf'. Another thruster parameter, 
accelerator current, is used to indicate when thi'^ conflition exists, since 
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it is primarily a functicr of propellant utilization. An override is, therefore, 
included in the loop to turn the main vaporizer power off at an accelerator 
current of 10 mA and to restore loop operation when the accelerator current 
drops below 6 mA. The override levels were Ctiosen to correspond to the 
experimental work presented in Ref. 11, 

B. DISCHARGE VOLTAGE -- CAIHODE VAPORIZE < LOOP 

This loop is shown in Fig. 5b. In this loop the discharge voltage 
is maintained at a level where minimum cathode erosion will occur. The 
electrical impedance of the discharge chamh r is a function of both the 
discharge current and tiie effective mass flow rate (M^). Empirical rela- 
tionships for these chamber characteristics are presented L. Figs. 11 and 
IZ. The experimentally determined value of K that was used for the beam- 
loop interaction was 0. 133 and was determined by using flow rate data. 

■ r ^ recent data on the 703 thruster indicate that K is closer to U. 1 Zb and 
. .-.y be a function of beam and discharge currents. The requirement for 
this K factor may best be seen in Fig. 13, where, for a constant ratio of 
discharge current to beam (5), the cathode flow rate decreases as the 
discharge and beam currents are increased. To keep the discharge v'oltage 
constant, one would expect the cathode flow rate to increase as the current 
is increased. However, the coupling from the main v aporize r causes the 
discharge control loop to lower the cathode flow' rate. The factor K may 
be determined with the following procedure. With the beam off, the cathode 
and main flow’ rates required to keep the discharge \-oltage constant over 
the range of discharge current required are measured. The sum of these 

flo-'v rates is defined as \1 . Then, with the beam on, the flow rates are 

c 

measured. Th( -. , for corresponding discharc* currents, K “ (M - m l/rhm, 

r c 

where rh and mm arc the tlow rates measured wit.i ihe beam on. 
c 

A flischarge supply provides the p ropell a nt - ionizing pow»T to the dis- 
charge chamber. As suggested in Uf'f. 1 (J, the discharge supply \' - J 
characteristics may be varied in order to stufiv the effects of power supply 
and thruster inte raclions. Several cha racte ri st ic s were invest ig at t?d , 
including constant current, \ to I r.itios of both i-) and I-) ti V 'A, and 
constant power. The effects of changine the V-I i ha raef e r i st i c s will be 
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examined in this report. The operating level of the controlled volta^^e sets 
the discharge po'ver for a corresponding beam-current level, usually at a 
ratio of 5:1, discharge current to beam current. This loop, therefore, 
controls the thruster efficiency. 

The discharge voltage is sensed at the output of the discharge supply 
and is compared with a normally fixed reference level. Voltage signals are 
acted on by a controller which adjusts a power supply to heat the cathode 
vaporizer. The characteristics of this controller are presented in Fig. 14. 
Thermal effects are included in the same manner as for the beam loop. The 
meicury flow again encounters a transport lag due to flow restrictions in 
the isolator, propellant line, and cathode. 

C. NEUTRALIZER KEEPER VOLTAGE -- NEUTRALIZER VAPORIZER 

LOOP 

This loop is shown in Fig. 5c. The mass flow rate through the neutral- 
izer is controlled by maintaining the keeper voltage constant. The flow rate 
is affected by the current drawn to the keeper electrode. The neutralizer 
characteristics used in this loop are shov/n in Fig. 15. The potential of the 
exhaust plasma is also indirectly controlled by this loop. It is a function 
of both the mass flow rate and the keeper current. In maintaining the flow 
rate, the keeper voltage is compared against a normally fixed reference 
level. Error signals, as in the case of the above loops, are used to adjust 
the mercury flow rate from the vaporizer. This vaporizer is also affected 
by the temperature of the on’/ironment, TE^. Very small delay times may 
exist due <-o the neutralizer isolator and cathode. 

In the laboratory, there is some interaction between the beam-current 
level and the neutralizer keeper voltage. There is little data available, 
and since the coupling is thought to be a function of the keeper supplv ''/^-I 
characteristic, this effect is not included in the model. 


reproducibility of Til . 
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IV. CONTROL LOOP ANALYSIS 


The computer model was used to perform some preliminary analysis 
of the thruster operation under both steady-state and transient conditions. 

In addition, the runaway conditions, under which a thruster escapes from 
the control loops, and the restoring mechanism were examined. Each of 
these three areas is described below. Experimental results of the thruster 
control-loop operation were used both to compare and to guide the modifi- 
cations to the model. Operation of the neutralizer was not included in the 
results described below. 

A. STEADY-STATE OPERATION 

Thruster operation has been observed to contain small oscillations 
(approximately 3 to 5%) in the beam current and discharge voltage. Some 
attempts at minimizing this type of oscillation have been presented in 
Ref. 10. The control-loops were, therefore, parametrically examined to 
determine conditions under which these oscillations may exist. The results 
indicate that if there is no delay time, there would be no oscillations over 
the range of loop gains used in the JPL laboratory. Introduction of delay 
times in both loops will cause small oscillations in the discharge and beam 
loops. The amplitude of the oscillations was found, w’ithin a small tempera- 
ture range, to be a function of the delay time. This analysis was performed 
for a low gain (1A/V) in the discharge loop and a high gain (110 A/A) in the 
beam loop. These settings are representative of current thruster operating 
conditions. 

An investigation of the effects of increasing the discharge-loop gain 
showed that the amplitude of the discharge oscillations increased rapidly 
as the discharge-loop gain was increased to greater than 2.0. These oscil- 
lations were coupled into the beam current. The investigation of variations 
in the beam-loop gain was not completed. Lowering the gain below approx’ - 
mately 90 A/A resulted in poor beam-current regulation. 

Identical delay times were used in each of the two thruster control 
loops. The beam-current oscillation waveforms are shown in Fig. 16 as a 
function of time for main and cathode loop delay times of 1. 0 s. The effects 
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of different le* els of temperature environment are also included. (In the 
figures, temj eratures are in millivolts output from a Type E (Chromel - 
Constantan) thermocouple). The amplitude of these oscillations are plotted 
against the temperature environment of the main and cathode vaporizer in 
Figs. 17-19. Nominal operating temperatures of the main and cathode 
vaporizers are also indicated in the figures. It can be seen that the oscil- 
lations are zero for a high-temperature environment (234°C or higher) and 
increase with decreasing temperature. A region of maximum amplitude was 
found to exist in each case with the amplitude increasing wdth larger delay 
times. The delay tim.es for the 700 series thruster were experimentally 
found to be about 1 s for the main and discharge loops. These delay times 
were measured (one at a time) by opening the loops and stepping up the 
vaporizer currents. The vaporizer temperatures and the controlled param- 
eters were recorded, and the delay time was taken as the time interval 
between the start of the increase in vaporizer temperature and the start of 
the change in the controlled parameter. 

This procedure is inexact, since there is an unknown time lag between 
the change in vaporizer temperature and the change in the thermocouple 
output. More work is required in this area, because the vaporizer flow rate 
is an exponential function of temperature and the actual temperature does 
start changing before the thermocouple output does. Thus, the delay time 
used is too long. This is pessimistic, since in closed-loop systems an 
increase in the time delay increases the tendency to^;^ard instability (a 
Nyquist diagram is rotated clockwise, or, in the frequency domain, the 
phase lag is increased proportional to frequency (Ref. 12)). 

The characteristics of the discharge power supply were also investi- 
gated to determine the effect of this supply on the oscillations. Three dis- 
charge supply characteristics are presented. These are: (1) a constant 

discharge voltage to discharge current ratio of 6 V/A, representing a linear 
approximation of a constant power discharge supply; (2) a constant current 
supply, and (3) a constant power output supply. The oscillations in the dis- 
charge voltage are compared in Fig. 20 foi constant current and linear dis- 
charge supplies. A 50% reduction in amplitude was found to exist in the 
oscillations whei. using the linear supply. Figures 21 and 22 provide a 
comparison of the oscillations using linear and constant po'ver output 
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supplies, respectively,. Variations in the cathode loop delay times had little 
effect on the waveshape. 

As expected, the constant power output supply minimizes the oscilla- 
tions. Figure 23a is a replot of Fig. 11 with the discharge supply character- 
istics added. For the same variation in mass flow rate, the discharge power 
change is less for the linear supply than for the constant current supply, 
and is zero for the constant power supply. Figure 23b is a portion of 
Fig. 6 showing these discharge variations coupled into the beam current. 

The flat slope of the constant power curve minimizes the coupling into the 
beam current. 

These results agree with experience in the laboratory. The supply 
used was ’’tuned*’ for minimum ripple, and when its V-I characteristic was 
measured, it was constant power over the voltage range of interest. 

B. THROTTLING 

The control-loop operation was examined during the transient condition 
that occurs in response to a change in thrust level. The response of some 
parameters in the model and some actual parameters is shown in Fig. 24. 

The differences between the simulations and the actual tests are due to 
(1) the difference in the magnitudes of the steps, and (2) differences between 
the model and the actual power supply available in the laboratory. The 
differences in the magnitude of the steps is roughly 2:1. Both the model 
and the laboratory results are higher than would be used in flight, since such 
steps would impose severe transients upon a flight power subsystem. The 
magnitude of the steps in the model was an exploratory attempt at under- 
standing the model, while the laboratory results when trying to duplicate 
the model were a practical limitation, since the thruster discharge would 
cease conducting if steps greater than those presented were tried. The 
w'ork on reconciling the differences is the reason for the author’ s opinions. 

In Figs. 24a and 24c the results of a step increase in beam and 
discharge currents is shown. The striking differences are in the magnitudes 
of the transient in the discharge voltage and the time required for the beam 
current to reach .he new value. Addressing the voltage change first, 
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^AVDM _-6.7 _ ^ ^ AVDA _ 1 _ , 

AIDM ■ 2.^5 " . AIDA " 1 " ^ 

,* 

Assuming a relatively constant DC plasma resistance due to time lags 
in the mass-flow rates the steps are in reasonable agreement, if the finite 
response time of the actual supply is considered. The time difference is 
thought to be caused by a hardware problem. The PP used for the laboratory 
tests is a modified 20-cm PP (Ref, 10). The major changes that had been 
made were rebuilding the discharge, screen, and accelerator supplies. 

The other supplies were left mainly unchanged, since the major effort at 
JPL concerned interfacing the TSS with a spacecraft, and plans were in 
process to procure new PPs designed for this thruster. Thus, the vaporizer 
supplies were power limited at a level lower than used in the model. Since the 
vaporizer power-temperature and temperature-flow rate characteristics 
are both exponential fimctions, the limitations in the power supply could 
explain the time difference. 

The difference in time response in the throttling down example is 
thought to be due to the difference in magnitude of the change in temperature 
required. The upward pulse at the beginning of the actual downward step is 
unexplained. During the testing, it was there most of the time. It may or 
may not be in the model; the computation interval used was actually longer 
than the length of this pulse, so the model may have missed this pulse. 

There is one other difference between the model and actual hardware. 
The step functions in the model are instantaneous, while, due to filtering 
delays, etc. , the actual discharge supply took a finite time to step the 
current. The model should be changed to be more representative. 

As the discharge-current to beam-current ratio is reduced (typically 
to 5:1), runaway occurs more often. The reason for this can be seen in 
Fig. 25, where, the higher the ratio, the farther the operating region is 
from the unstable region. The dividing line between the stable region 


AVDM = model change in discharge voltage 
AVDA = actual change in discharge current, etc. 
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and the unstable region is the locus of points in Fig. 6 where *he slopes of 
the constant discharge power curves become negative. Note in Fig. 25 that 
a constant current discharge supply has the shortest path to travel to the 
unstable region in rase of a perturbation in the discharge. It would move 
vertically, while other V-I characteristics would move parabolically or 
with a smaller slope. 

C. RUNAWAY CONDITIONS 

If thruster operation is shifted to the portion of the curves in Fig. 6 
where the slope is negative, then the loops will no longer function correctly 
and the main vaporizer power will be supplied at the maximum level of the 
supply. This condition has been referred to as "runaway" or "low-mode" 
operation by various investigators. Several thruster parameters can be 
used to indicate when this condition occurs. One of the most sensitive 
appears to be the accelerator current. A controller that operates as a func- 
tion of the accelerator current was therefore introduced in the loops 
(Fig. 5a) to restore proper control-loop operation if runaway occurs. This 
controller turns off the main vaporizer when the accelerator current is 
greater than 10 ^!A and re res loop operation when the accelerator current 
falls below 6 MA. These current amplitudes were suggested in Ref. 12. 

In the model, using a constant current discharge supply and a 1-s 
time delay, if a momentary runaway is forced it will not restabilize as is 
shown in Figs. 26-29. Figure 2b shows the beam current and discharge 
voltage variations as a function of time after the initial runaway. This run 
was terminated by a computer run-time limitation without ever stabilizing. 
Figure 27 is a variation of a phase plane plot showing beam current and 
main vaporizer flow rate as a function of time. The returns from the high 
flov rate conditions were caused by the runaway trip switch. Figures 28 
and 2^ are similar plots of these parameters obtained under the same con- 
ditions except that the delay times were set to zero. After some large 
oscillations, the run did stabilize with small oscillations in both beam 
current and discharge voltage. 

Figure 30 shows the computer outputs when a discharge supply with a 
6 V/\ slope and delay times of one second were used. The beam current and 
discharge voltage stabilize rapidly, followed by oscillations of less than 1%. 
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V. CONCLUSIONS 


A hybrid model of 30-cm ion thruster and treinsistorized power 
processor was formulated and programmed on a digital computer. This 
model was exercised over a range of parameters and the result led to the 
following conclusions; 

(1) The control loops, as presented, are stable for normal values of 
loop gains when propellant line delay times are zero. 

(2) Small oscillations will exist in the loops when representative 
thruster delay times are included in the mcdel. The amplitudes 
are on the order of 3% as normally present during thruster 
operation. They were found to be a function of delay time, loop 
gains, and the thermal environment near the vaporizer. 

(3) Discharge supply characteristics can affect the control-loop 
oscillations. A constant power discharge supply is most 
effective in minimizinv these oscillations. 

(4) The model quantitatively predicts transient operation. 

(5) Recovery from runaway can be aggravated by delay time, 
and to a greater extent, discharge supply characteristics. 

A constant power discharge supply was found to improve 
recovery characteristics. 

The results of this initial work appear to be a useful tool in understanding 
control-loop phenomena. Continuation of this type of effort should provide 
an optimization of the runaway- recovery techniques, and eventually provide 
a model able to duplicate parameter variations that may occur during a 
mission using solar electric propulsion. 
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Table 1, Normal thruster operating and measured parameters 


Parameter 

Operating point 

Discharge voltage reference 

37 V 

Discharge current reference 

5 A 

Beam current reference 

1 A 

Delay time for main vaporizer 

1 s 

Delay time for cathode vaporizer 

1 s 

Cathode vaporizer heater gain 

1 A/V 

Main vaporizer heater gain 

110 A/A 

Magnetic baffle 

40. 5 ampere-t irns 
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Table 2. Ion thruster description 


Thruster 

Baseline 

Grids 

Electron 

Baffle 

Neutral- 

izer 

S/N 404 

Ref. 9 

Replaced flat grids 
with curved grids 
0. 5% compensa- 
tion 0. 1 5 cm dia 
accel holes 

Increased 
to 5. 87 cm 
dia 

No 

change 

S/N 703 

Ref. 10 

0. 4% compensation 
in place of 0. 5% 
compensation 

Increased 
to 5. 40 cm 
dia 

400 
series 
with 
electri- 
cal iso- 
lator 
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1. 700 series thruster 


\2 kAOLAL '^RMANENT 12 AXIAL PERMANENT MAGNETS 



Fig. 2. Discharge chamber schematic 
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Power processor grounding 

































Fig. 5a. Block Diagram - main vaporizer loop 



Tfl 


Fig. 5b. Block diagram - cathode vaporizer loop 



T13 

Fig. 5c. Block diagram - neutralizer vaporizer loop 
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ORIGINAL PAGE IS POOR 
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PffOPELLANT FLOV' RATE, g.1> 

Fig. 6. Beam characteristics 



Fig. 7. Main vaporizer heater supply 
characteristics 
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Fig, 8, Vaporizer time constant 
(bell jar tests) 
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Fig. 9. Mass flow rate vs temperature 
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Fig. 14. Cathode loop heater 
characteristics 



Fig. 15. Neutralizer keeper 
characteristics 
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Fig. 16. Typical computer runs showing the effect 
of the environmental temperature 
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ABEAM CURRENT, mA 



A DISCHARGE VOLTAGE 




8 10 12 U 16 18 20 

TEMPERATl RE, mV 


Fig, 19. L.oop oscillations vs temperature 
at typical beam current (2.0 A) 
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A BEAM CURRENT, WA 






0 2 4 6 8 10 12 


TIME, » 

Fig, 20. Comparison of discharge 
voltage ripple 
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time, $ 

Fig, 21. Oscillations using a constant 
slope supply 
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Fig. 22. Loop oscillations using a 
constant power supply 
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Fig. 23. Discharge supply characteristics 
plotted on thruster characteristics 
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(c) 



Fig. 24. Comparison of simul 
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and actual step responses 








REPRODUCIBILITY OF THE 

ORH’INAL page is poor 



Fig. 25. Discharge characteristics 
showing the effect of 
Id/^b ^^tio 



Fig. 26. Runaway recovery characteristics 
when using a constant current 
discharge supply 
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APPENDIX A 

COMPUTER -PROGRAM LISTING FOR THE 
ION THRUSTER MODEL 
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REI^RODUClEU-riY OF THF 
ORIGLNAL PAGE IS FO)H 

-AbGfA Ci5L*TKAN 
-XuTiO CbSL* TKAi\*CSSL 
PkuokANi lui\ 

CONiP.ENT 700 StMEo TrlRuSTtR 
COi'iMcM' TnIG PkuGKhi'I PuKt'-S rli^Au SnAPE 
COi'inCiM OuAi^ TiPtC IStL tOk MAi fs i_OoP 
^OP.Hc.M ucMu» TiMt. 1 :>tC rCK CATnOOc. lCOP 

LUi*i>iwi\T ouoAl. PC/»k£K ;jOPruY 

Cur-iMtM fli*i£ CUi<i)T 

CCi\o1Ai>iI 1 A\>i = £ 10 • o » rAu2- j 0 *0 * TAo3-21C«0 

CuMMti\T I u.^PcXA I jkc.-HuOT 

CONST Ai\l AuPriAi = 0 • 33 I o 1 = -7 •0533 

constant ALPnA2 = 0*35#i2=^-3*7625 

Constant alPha3=o« 33* ci3=-7*o693 

COl'iMcNT i-iLATtA Nc.S 1ST ANC£ = C0NST OnX 

CONSTAi'^1 A = £*6S 

Cui'tpic.uT OibCnAAuw CnAMoEK Cdmi^ACTEk I S f I CS 

CONSfAi'iT vAt\Ck=^7«o 

CuHiMLi^T y.ASi' OF Y 

CuN^j r Ai'i T .'’ih’j- 3 • j 3 

CoN.Mci\T Dt-Mf'i CJ^M^c.NT 

CONSTAi^i ioKErO-=i*3 

Cui'iMLNT GaIim 

C'-'NS T AN T GA i N ir\ = 1 • v. * G A I n 2 = ~ilC’»0» uA Ii't3K=i*0 
CGNyiENT lONTKOL pOINT 

CO(\i T AN T ^-ii“c»3» v-i3 — 0*3* C4.1“<'*w» C^3 — 4«0 

(_Oi If-'C-N 7 ivcu i fxAL i 4Ct\ KC-EKLrt CukkcNT 
Constant ii\si-<=o»o 

i N i- >j T k 'll- 1 1. i„i\ s. — v0i-Tm\jl 
CO i'nS TmiS T V I »N.f\ — i • - 
CO(‘i*itNT LuV 1 i\0.^'''i£N t I C-iSP c.\ A T jRc. 

C-ONSTAi'mT 

CON:>TAi\l 7c^i=i3*o 

CONSTAinT TcJ=i3«0 

COMiXci'tT Vi>-IiN Flow a f T c.-'iOA T i uN 

C0.'.S1A(\T \ = u.i33 

CONSIANT oi=l#o -19 

CoNS^«n1 Yuu=-4»3 

constant Y01=5*3 

OmlL y.hiN 

SLcPtNi iVAlN 

INi T lAL 

pAot Eject 

Fi ILL I »V‘Ai\C»lo* iMLCEi_»V''t\»YLT'^J 
II I 0 L oPiC » oPf'i » vj i*ir. tIl»IE»Ij 
lilLE TC.'4.1 »Tv. 21» F L AkC » ki T i M C » I AkC 
N 1 = 1 • 0 

COiMMlNT ASS*oN CA]m vAr' ntATEk INH'Al CukRENT 
iilC=i»^ i iii=i«2 

COViMlnT ASSIuN Vmt ncAicr^ INIiiAL CUKRENT 

I^IC~0»c i lZi~0«to 

Coi*iyLi\T mjSIoN NlvjT V AF ri:.ATt.k IimITiAL v-okkcNl 
IjIC=w»o j> i3i = 0»J 
1 AkCk 0»= 7 , b 


36 


JPL Technical Memorandum 33-755 



GAlN3 = bAli\3K 

ACCN-O.U 

tNL) 

DYNANiiC 

VAKlABLt ;=O.0 

CONiMENT .\C' CHA'^GE Ih AkC CUkRENT Ii\ TciIS CASE 
AiiRAkC=0 • o 
1 AKC«n= 1 AkCKvj-AciKAKC 
1 bKtP = i DKtr 0 

CU/iNiEM bAlnuOt. Vht uMiN CnAfWGc 
AbKo i = S TEk I • 0 » I ) 

AbK(j2 = bT tKlHUU^Otl) 

Ab«Gj = STEH(ouvJ»J*l ) 

<'akoa = ST tP I dCO* C * T ) 

Ab RGb = S1E^'( 1000«G»r ) 

GA i i\ i = «jM i N 1 k-f AoRGi + AbRvjZ + AuRo3 • • • 

AbKOA +AbRG5 

CUM^-EM calculate c-ACri VMp nbATER CURRENT 
MACkG i'>iA'^KO mEATER Q 
yACRu RtLAoc-L LltLi: 

I F ( u 1 A ) •Gc • 2 • 0 ) oO TO LI 
Q(l)=0(b) j> u(^)=0(6) 4> u ( 3 ) =-0 1 7 j 
GU TO L^ 

Ll»*<J(i)=0(o) i •J(t)=^(‘y) 2> 013) =0(10) 

L2 • • lON f 1 Nut 
MACRu tND 

riEATER iifi^»I3*('ti*IliC»I2iC»I31C»»*» 

.V 1 =K'>' i i i rt2=R*i2*'^2 X 'rt3=R-* I 3*‘''2 

STlJi = TL»A(Ai} i o7L*t = TU>N(ir»2J i oTD3 = TDft'(w5) 

lF(iM«oc,^.u) GO io l3 

TOAi^bfUj. X. TD/ii=ofL-^ X TDA3=Gi03 

GO Tu L*+ 

L3 • • 1 Ox i>.'= TOX 1 X TuX2G = T0a£. x T0X30=7Dx,‘> 

Tui = olJA i> T0^ = -T04 i T0-> = iT0? 

lmlL iiNil luiXlufTOxtOfiOXjJU’^Oi'TOZjTOG*! Oa i*T0x4.» 10X3) 
L A • • CuN I I NUc. 

I L,W Y ioi'< 

Ot K 1 ViA T i Vc. uui< 

SuOi''i=INTto{ ) 

CIMEkvAL Ci=C.2 
tNU 

LUi'..'*iCts I \/ At'ur. 1 ^Lrl r C-Mr LKA I oKL 
I Cxi = ruxi+ I tx 
i C^; i = ' 0X2-^ T E^ 

Tc3i-=TuX3-^Ic3 

Comment ma^x Fww.v i-^mTu 

SMC = ExP (Bi ) kMLi'nAi* T*. ii < 

SMiX = c.a(- i o2 ) »c.Xr^ ( Al»- nA2* 1 C2 i ) 

SMN = EXP (5jJ ’'Exf (Ai_rriAJ^TC5i; 

LOMMEhI OtAO ( I .t A,\ T A i i 0i\ «_*vu j 

1K(N1»ol«2«oj oO To l>u 1 
Yii = Si'''C X Yii = Gi*lC T r xi-oMC 
Y4l=SMC i Y^i=oX^ 
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Y6i = SMC S Y71=SHC i Yol = 5.viC 

Y91 = SMC i YiUl = byiC 

Yi2 = SMy i> Y22 = 5i'ii'i i> Y32 = oMM 

Y42 = bMi*i S Yb2 = iy.vi 

YlJ = by/Y i Y23 = bWt'4 i YJJ = bMiN 

Y43 = bMi\ S Yb3 = bnu 

Yti^byC b Yt2 = Sni*i i Yc3 = ;3yi'l 

bSMC = Y£l i boMh-^YEL $ bbXN = YE3 

GU TO 0D2 

l)L)1* *^uNT li'^UE 

MACko NiMCkU UtLl U 

MACRO kEDtFikt Ym i . Y 51 * YC 1 . YDl * Y L 1 
MACRO KlDLFINl t i- 1 , YGl » Yhi ♦ Y i 1 » Y J 1 
YAl=0(i) » Ydl=o(2) b YCi=U(3) 

YDi=U(A) i YEi=u(p) 

YFl=Oio) * YGi=u(7) i YHi*0(3J 
Yli=J(v) * Yji=Uliu) 

U 1 11 ) = Y Ei 

U(i)=U(12) * U(<i)=TAi i u(3)=Ydx 

U( A) =YCi b 0(5) =YOi 

U(6)=YE1 b 0(7)=Y)-i b U(d)=YGl 

0(9) =Yrti b 0( IJ ) =Y1 1 

MACRO END 

MACRO i>'iACRO DEL2 00 

MAv.RO Kc.Dt.rlNE YA2 »Yo2*tC2»yD 2»YE2 
YA2 = 00(i) i- YD2=uOlii) i YC<: = U0(3) 

YD2=00(A) i YE^=Ou(3) 

UU ( 6 ) =YL2 

U0(i)=ou(7) * Oo( 2)=YA2 b Ou(3)=t62 
UU( A)=YC2 S 00 ( 5)=Y02 
MACRO END 

MACRO MAlkO 0EL3 000 

MACko Rt.Dt.rlNE Y A3 , Y lO ♦ Y C 3 ♦ Y03 , Y E 3 

YA3=000(1) b YS3=OUO(2) i YC3=JOo(3) 

YD3=OOu( 4) b YE3 =uOO(5) 

000(6 ) =Y£3 

00c, ( 1 ) =000 ( 7 ) i 000(2)=YA3 k OJJ(2) = Yb5 
0001A)=YC3 b 00o(5)=Y03 
MACRO LND 

DElI Yii»Yii»Y31»Y4i,rPl»Y6l»r71 tYci *Y91»Yl01»o SMC t Oi'iC 

DEL2 Y12»Y22»Y32»YA2»Y52» boMM » bM -i 

DEL3 Yl.J»Y2b*Y33»Y‘+3,Y53» SSMiN »SMN 

DD2. .CONT iNoE 

MC = SSMCi-K*-bbMM 

MT=SSMC+SSMM 

,MN = SSMi\ 

CCMMtNT TABLE OE VAPORIZER TcMPlRAToRE ChANoE CORRobFONo i No 10 
TABLE T DW » 1 ♦ 6 * ••• 

0*0(1«5« 3«5y 6«0i 9*u» lG«6t ••• 

C>Cy2*Ot 5*5* 10«5*lf*Dt23*0 
COMMENT AKc CmAimolR CHAR Ac TER I 5 T i Cb 
TABLl Ar<CMC » 2 * 7 t / » ••• 

!• 7d*l*'7t^«Js«^*At2*b>2« 7Dy3*5t • • • 

3iO»3*u« 7»vJ» 9«o» ll*Ofl3*0»13*0»*** 
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61*3»37*C»^6«3»20«u»lb*C»10aO>6*f>»*a« 

76*0>32«'Jt3^*Gt2o« l •• 

90*G*67«Ct^daOy37«G»3G*3*23*3y21a0t*** 
l02a‘Jy7/«Gy3t5a3y^3aC'y3/»Cy3Ga7 2<a0yaaa 
lilas^yOOaUyt)7 •'jy33a 7 f^3a3y37«'Jy 32a0y a a “ 

116aQy93aO» 7^aC't60aGy?Ua3yH<.a3y 3 7 «0 
COMME.NT bcAH CHAkAC T tk I :j i i Co 
TAbLb iuMry2y22y7yaaa 

0aGyAa0ybaUybaGy^a*jytiaCy^aCyi0a0ylia0yl2a0yaaa 
13aOylAa0yl3a0yl6a'^yi I aOylOaOyl 'a0y20a0y^la0y^^a0yaa « 

23aCy2Ha0ya aa 

lllaOylHbaUyid3a(jy<^22a0y^b'^i 0f29b«0y 333a0yaaa 
GaUyUa3 7 yv>‘aH9yUa6C'yCab6f‘Ja 7ij0 7‘+»C« 75*0a 7^»Oa7 iyaaa 
0a66y0ab0y0a32yba^9y0a23)0ajby 'Oal3y“Ga‘*7 y—J a07 y'~'Oa93y«aa 

la23y~labCyaaa 

CaL>yUaA2yu«b<yyUa&3yOa7^yUa6<.iOa6byUa 92y0a laOO* a aa 

Ca9'yyC'aVbyOa>iyUabA»’jab6fO» >ly0a3ba0al3y“0al4y — 0a3byaaa 

0abG(~VJaOb ya aa 

’•JaL'ytJaH2yUa:>byUa6/yOa6O''*'>9iyla00yla09ylai3yialbyaaa 

la,^Uyiaioyiai3fXaU4»'J.-7<.'‘»9ab4yGa73yUa64yCab0y0a33yaaa 

£. 3-9\J a 1 b y a a a 

OaCyUa^byL/ab / y0a70y0aO^t0a'^7yla0Oyla 10yla^byia33yaaa 
la3 / ylaHl»la‘t2fla‘+i»ia09yia33yia^.6yi a ibyiaObyOaVbaaa a 
OaObyUa 74yaaa 

0aOyUa^7 yUabl tOa 7*+jL'abb»Ga9 7>ia I2*la^Ayia3*y* laH^yaaa 
laPlyiabb*lab0ylabiyiab3yiaOlylabb»l abbyla50»la43yaa a 

labbyla^byaaa 

JabyCa4^ yOaOl yOa 74y0aOby'Ja'7 / ylai2«la2byla36y la49yaaa 
laOC y iabV* 1 a 7o yiaOJy la :>Cy ia 7by2a0')y^l a'j4a2a'j6y2aC9y a a a 

2a lUy^av>9ya aa 

OaU*Ga47yUablaOa/^aOaOba^a'7/yial3yia<^7alaHlalab*4aaaa 

iabb»ia77yia90yia^'7yi:aljy2ax:Uy^a30y2aH0y<.abi92abiyaa, 

2a71y2a 

COyiMc-iV 1 NL^T <LLi CK ChAKhCIlK 1 ijT 1 cs 
TAdLL iNKMN*2a6y‘+a aaa 

OabyC'a/aQaOy 0a7» laQy laiy aaa 

Qa31yGa^OyCa4lt0a4Oy aaa 
13abyx^abyl3abyljabyl.?abyl3a5y aaa 

13aVJal3aCyl3aCai3aOal3aQai3a0y aaa 

12abyjLCabfi2abyi^a'-^tl4a^yl^aby aa a 

12a>^yi2aOyx2a0yl2aUyl2auyl^aU 
NtO= 1 a U 

IF(l\t^aLTaOaO) oO TO LNLo 

CO^^l•^LM NtUT Pu'aLK buPHLY CHAPC T lH I S T i Cb 
YOI AKC=7aU 
SLuPMC = b a ^ 

I AbC - I Yv^lAiXC aOaUOi y20 yFLAXL y iAKCk'*‘a a a 

( V AkCK-AKC‘*llw I i»1L ylAKv-) )/bLUKMLyOaOi) 

LN Lb a a C OIM T I i\UE 

VAKC=Ai<CMC ( MC a 1 AaC ) i AkC = VAKc * i AkC 
I&=IoMT IMT aPAkC ) 

Y01M; = Uab 

SLuPl''^^ = 3uaU 
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lNK>=iHKL(YOiiMK.»0*(Joi*iU»>’ LNl^» lNKr^ + ** • 

( VKK.n~lNKMr<i ( i » i^it\ ) ) / SLor'i'iK. »0 • w i ) 

VNN. = i\fvMiM ( Ii\kfi*!:\) 

COWMtM CATnUDL VAr'OKlz.t.K CUiN rixOLi_c.r< 

Vl=(Cii~C^l)/l2 1 ) + VAkCK 

V2 = (C21-Cli)/(2 • 0*oM ii\ i ) ■•■ V ArsCf< 

PKuCLjUKAl 1 iii = VifV2»VAK(_*VAr<CK»'wil*C^:l»0Mli\l ) 
IF( VAkC.lt. V lJ Ili=Cil 
1 f- ( ( VAKC.Gt.Vl ) . Ai\l3. (vAkC.Lc.V 2)) ill = ... 
GAiM»VARC-oAINl»VA>^Cr<-r (C11+C21 ) /2.0 
lF(VAKC.GT.v2) lii=C2i 
c.Nl' 

COi'^MLi^I N'AIN VAPORIz-hK CGNTR0»-L lX 
01 o= io-ioKcF 
DIb2=lv_i2“C^^ )/oMii'4t 

Pi-<OCc.UUr<AL( Iii=u»iD^*L‘l o»Ci2*C22»uAli>t<:) 

IF( (Jlc.ou.U.O) .MNu'.(k^lb.Ll.0Io2) ) 121 — ... 

lA I N<:*UI D+C22 
1 F ( 0 1 a . LT . G . C ) i<ii = C^2 
1 F ( u I b . oE . J i o2 ) I2i=Cl2 
lNL' 

COMi'lt-Kl NduTn.ALiZLrt V APUi-^i 2 Crv. CU T R OLL LK 
VM^l = lCi3-C<i3) / ■•■V.'^NK 

VNA2=(C23“CiJ) / -r V i>*nK 

P KOC ULUKALi ijl = VKS.»vA\l ♦Vi'ts^.»Ci3»C^^*vjAif«a»vNN,K) 
i f- I vnA . LT . Vi\k 1 ) i3i=«^ia 
Ir( lV^^^.uc..VNAi).M(^lu•(\/l^l^.LT,Vl'i^2) ) i3l = ... 

GAIN3*VN«-GA1 Na * VNnK+ lCl3 + C^3)/2.C 

I F ( ViNi^.ul . VfsiKt ) i3i-CziJ 
C.NU 


COi'^MEiNT ACClL CukK£.i'*T I M lim^uPT 1 0i\ 

lU = ^l*y.T / iMhli*- 1 . ul- 3/ 3oOC. G * YLTAu=io/Iu 

TAbLt MCYc.lA»i*i'j»..* 

C.l*C.2»G.-3*0.4yC.D»0.btO. / >0.b)0.v*x.Gv... 

0.L27*a.019*'^t0.bl3y0.011^ tc.G0'3b *... 
C.00b6»L.uC3»0.CC37,0.C0^V*C.Cj22 
lACCtL = ACY£:TA( t'£lAu)*iD 
1 F ( 1 Av^C C.L . oc. . 0 . 0 1 0 ) oO To A„li 

IF ( lAO^_LL.Lt_.0.LVi^fc>.'Ai\i^. i<^OoC-L.i-T .L.O^C.^'i'ic/.Av.v.i'I.Lo.l.O ) oG lo 
ACi i 

lF(iALLt-L.'jL.3.00o.Mi\G.lMoCcL.LT .<.>.GxO.'->i\L).ACv_iN.t'Ji.U.0) oo TO 
AC I 3 

I F ( 1 ACCLL . L T • U • OCb ) ou TC AC I 2 
AC 11. .121=0.0 3) ACC.\=i.O 
GC TU AC 13 
«C12. .ACCN=O.0 
ACi 3 . .Cunt INUE 
I F ( T . L T . 2<J . u ) Gu Tu la2 

IF(T.Gt.io0.0.ni\u.T.i_T.200.0) oO Tu LW2 
IF(T.Gl.3o0.0.hNu.T.lT.hG0.0) uu Tu l'a2 
IFIT.UL.doo.O.ANG.T.lT.oLO.O) Gu Tw i_a2 
IFlT.Gc./oo.O.ANUfcl.ul.oOj.O) UO Tu la 2 
IF(T.ut.^o0.U.ANG.>.u.r.i00u.C) oC) To LW^ 

I F ( T .ut . 1 io 0. 0. AND. T • L T . 12u0.0 ) ou To 


• • • 


• • # 
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'iDUM=l.G » GO 10 LWi 
LW2 . .CUNT iHUt: 

OUT TtVAKCf Ib»IACCt:L.vNK.»YtirAj,CI 
OU T SKC » SMM * 6MN 

OUT TCli*TC2l »FLAkC*MI .MC*IaiK» ^1 

PAot SK.iP»itLl 

LW 1 • .CUNT 1 Nut 

IF ( 16.Lt.U.0) GU TO LP 

IF ( T.Gt.AcO.C ) GU TO LP 

N1=2.U 

END 

terminal 

Lb#. CuNTlNut 

END 

END 

COMMENT FiKoT OkucK l«o TKMi>icjrLrv FuiS'^.TIGN 

SEGMLN T 1 NT I T l>A 10I-' » 1 DA^uH »1 DApO^' »TDiM »1 Lcr.TUJp «••• 

TDA1P»TDA^M»TDApP) 

initial 

N N = o • 0 

V = U.v. 

TOAiu=TUAiP » TuA4U=IuA^h * TuA 3 >j. = lu A 3P 
TDiQ = TDiP i TD^u=Tu<:P j> fLPu=[UP>-' 

TDA luu= TDX lOP J> T DAlO’---TDa^OP f Da p Ou- T JX pGk 
END 

DYNAMIC 

VARIABLE V = «-'.0 
IF(V.Gt.0«2) Gu TO L6 
DtKl VAT i Vt SINT 

TDAlu = KtALPL'. TAui» lu'l>*»TOAx''Ci.') 

TDA2u = f\cALPLi TAU2 » 1 u2u* TUA^UC) 

TDA3- = KlALPL( TmU3» Tu3^* »TDA30Q) 

CINTERVAL CIP=0.2 

NN=1.U 

END 

END 

terminal 

L6. .CONT INJE 

TDXlP = TOAlCi 5 > TDX*iP = TDX2u i TDa3P = TDA30 

TDlP = TDiO 3) TD^P = T‘J«iD j> IDpP = Td30 

TDX1UP= TDXlOO i TDX^Gk- TDA^CQ » I DA p OP = T DX 300 

end 

END 

END 

-FIN 
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